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Abstract. The properties of Dirac electrons in a magnetic superlattice (SL) on 
graphene consisting of very high and thin (5-function) barriers are investigated. We 
obtain the energy spectrum analytically and study the transmission through a finite 
number of barriers. The results are contrasted with those for electrons described by 
the Schrodinger equation. In addition, a collimation of an incident beam of electrons 
is obtained along the direction perpendicular to that of the SL. We also highlight the 
analogy with optical media in which the refractive index varies in space. 
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1. Introduction 

During the last five years single-layer graphene (a monolayer of carbon atoms) has 
become a very active field of research in nanophysics [H [2]. It is expected that this 
material will serve as a base for new electronic and opto-electric devices. The reason 
is that graphene's electronic properties are drastically different from those, say, of 
conventional semiconductors. Charge carriers in a wide single-layer graphene behave 
like "relativistic", chiral, and massless particles with a "light speed" equal to the Fermi 
velocity and possess a gapless, linear spectrum close to the K and K' points. One 
major consequence is the perfect transmission through arbitrarily high and wide barriers, 
referred to as Klein tunneling. 

One of the most challenging tasks is to learn how to control the electron behavior 
using electric fields in graphene. This task is made complicated precisely by the Klein 
tunneling according to which Dirac electrons in graphene can tunnel through arbitrarily 
wide and high electric barriers [3]. 

Alternatively, one can apply a magnetic field to control the electron motion. It 
was shown in numerous papers that an inhomogeneous magnetic field can confine 
standard electrons described by the Schrodinger equation [U El E]. The question 
then arises whether it can confine Dirac electrons in graphene. Up to now semi- 
infinite magnetic structures, that are homogeneous in one direction, were considered 
and made the task simpler by converting the problem into an one-dimensional (ID) one 
[3 El El HO, EH [El EH [HI [15]. In particular, a magnetic confinement of Dirac electrons 
in graphene was reported in structures involving one (7] or several magnetic barriers 
[SI E] as well as in superlattices, without magnetic field for some very special values of 
the parameters involved [IB] . In such structures standard electrons can remain close to 
the interface and move along so-called snake orbits [5] or in pure quantum mechanical 
unidirectional states [I]. 

Given the importance of graphene, it would be appropriate to study this magnetic 
confinement more systematically. We make such a study here by considering a magnetic 
Kronig-Penney (KP) model in graphene, i.e., a series of magnetic ^-function barriers 
that alternate in sign. This model can be realized experimentally in two different ways: 

1) One can deposit ferromagnetic strips on top of a graphene layer but in a way that 
there is no electrical contact between graphene and these strips. When one magnetizes 
the strips along the x direction, cf. Fig. [H(a), by, e.g., applying an in-plane magnetic 
field, the charge carriers in the graphene layer feel an inhomogeneous magnetic field 
profile. This profile can be well approximated [22] by 2B zoh/d(x 2 + z^) on one edge of 
the strip and by —2B z h/(x 2 + Zq) on the other, where zq is the distance between the 
2DEG and the strip, and d and h the width and height of the strip (see Fig. QJb)). The 
resulting magnetic field profile will be modeled by two magnetic 5 functions of height 
27iB h. Such ferromagnetic strips were deposited on top of a two-dimensional electron 
gas (2DEG) in a semiconductor heterostructure in Ref. [23J. 

2) It was recently shown that local strain in graphene induces an effective inhomogeneous 
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magnetic field [24J (Fig. Q](c)). When one puts the graphene layer on a periodically 
structured substrate the graphene at the edges of the substrate becomes strained and 
the situation can be described by a magnetic 5-function profile such as that shown in 



In a quantum mechanical treatment of the above two systems the vector potential 
A(:r) is the essential quantity and, within the Landau gauge, A(x) is nothing else than 
a periodic array of step functions. The Hamiltonian describing this system is periodic 
and consequently we expect the energy spectrum of the charge carriers in graphene to 
exhibit a band structure. The advantage of this magnetic Kronig-Penney (KP) model 
is mainly its analytical simplicity that provides some insight and allows a contrast with 
the same model for standard electrons j27j. To do that we adapt a method developed in 
optics, for a media with periodic in space refractive index. This optical method is clear 
and very well suited to the problem. Incidentally, there are many analogues of optical 
behavior in electronics, such as focusing [TTJ, [18l [19] , collimation or quasi-lD motion of 
electrons and photons [TBI EQl SI [HJ , and interference [21] in a 2DEG. 

The manuscript is organized as follows. In Sec. 2 we present the method and 
evaluate the spectrum and electron transmission through two antiparallel, ^-function 
magnetic barriers. In Sec. 3 we consider superlattices of such barriers and present 
numerical results for the energy spectrum. In Sec. 4 we consider a series of ^-function 
vector potentials and our concluding remarks are given in Sec. 5. 



Fig. m 




Figure 1. (a) Layout of the system: a ferromagnetic stripe on top of a bilayer graphene 
sheet separated by a thin oxide layer, (b) Magnetic field and corresponding vector 
potential at a distance zo = 0.1 under the stripe for parallel to it. (c) The graphene 
layer on top of a periodic structured surface. 



2. Characteristic matrix for Dirac electrons 



An electron in a single-layer graphene, in the presence of a perpendicular magnetic field 
B(x), which depends on x, is adequately described by the Hamiltonian 

H = v F cr- (p + eA(z)), (1) 
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where p is the momentum operator, vp the Fermi velocity, and A(x) the vector potential. 
To simplify the notation we introduce the dimensionless units: £b — [^/^Bq] 1 ^ 2 , B(x) — > 
B B(x), A(x) — > Bo^A(i), t — > Ub/vf, r — > £sr, v — > vfv, E — > -Eo-E 1 , — >■ 
E u(x), Eq = hvp/f-B- Here is the magnetic length and t the tunneling strength. In 
these units Eq. (CQ) takes the form 

H= ( d x - idy + A{X) \ 

\d x +id y -A{x) J • 

Then the equation H^(x,y) = E^(x,y) admits solutions of the form 

ipi(x,y) \ 
v *l>n{x,y) ) ' 
with i[)i(x,y),i/)ii(x,y) obeying the coupled equations 



*(s,y) 



(2) 



(3) 



— - i— + A(x) 
dx dy 

— A(x) 
dx dy 



ipn + E^ = 0, 



i>t + EtPjj = 0. 



(4) 
(5) 



Due to the translational invariance along the y direction we assume solutions of the 
form ty(x, y) = expik y y(U(x), V(x)) T , with the superscript T denoting the transpose of 
the row vector. For B(x) ~ 5(x) the corresponding vector potential is a step function 
A(x) ~ B(x). For A{x) = P constant, Eqs. (jlj) and (jSJ) take the form 



dx 



V = iEU, 



U = zEV, 



(6) 
(7) 



Equations ©-(Ej) correspond to those for an electromagnetic wave propagating through 
a medium in which the refractive index varies periodically. The two components of 
$?(x, y) correspond to those of the electric (or magnetic) field of the wave [28, 29J. 
Equations (jSJ) and Eq. (jTJ) can be readily decoupled by substitution. The result is 

d 2 Z 



where Z 



dx 2 
U,V. 



+ 



e 2 - (k y + py 



0. 



If E 2 



E' and {k y + Pf 



V eff , Eq. 



Schrodinger equation for a standard electron where V e ff(k. 



x 



(ky 



(8) 

reduces to a 
-f P) 2 can be 

considered as an effective potential. Taking #o as the angle of incidence, we have 
k x = -Ecos^q = [E 2 — k 2 } 1 / 2 and k y = Esm9 are the wave vector components outside 
the medium and k' x = Ecos9 = [E 2 — [k y + P) 2 ] 1 / 2 is the electron wave vector inside 
the medium and 9 = tan -1 (k y /k' x ) is the refraction angle. This renders Eq. ([8]) simpler 
with acceptable solutions for U and V 



U(x) = A cos (Ex cos 9) + B sin (Ex cos ( 



(9) 



V(x) = -i {B cos (9 + Ex cos 0) — A sin (9 + Ex cos 9)}. 



(10) 
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For future purposes, we write U and V as a linear combination of Ui, U 2 and Vi,V 2 : 



dVy 

dx 
dU x 



+ {k y + P)V 1 = iEU 1 , 



(ky + P)U l = iEV U 



dV2 
dx 

dU 2 



(ky + P)V 2 = lEU 2 . 



(ky + P)U 2 = lEV 2 . 



dx dx 
We now multiply the equations of the first row by U 2 and Ui, respectively, and those of 
the second by V 2 and V\. The resulting equations lead to 



5? = u[v 2 + UxVi - v;u 2 - v x u' 2 = o, 

ax 



where D 



(ietZ) and 



u 2 v 2 



(12) 



(13) 



Equation ffT2l shows that the determinant of the matrix f|T3|) associated with any two 
arbitrary solutions of Eq. (|SJ) is a constant, i.e, D is an invariant of the system of 
Eqs. (jTTj) . This also follows from the well-known property of the Wronskian of second- 
order differential equations. For our purposes the most convenient choice of particular 
solutions is 

U 2 = F(x), 
V 2 = G(x), 



Ui = f(x), 
Vi = g(x), 



(14) 



such that 



f(0) = g(0) = 0, F(0) = G(0) = 1. 
Then the solution with U(0) = Uo, V(0) = Vo, can be expressed as 

U = FU + fV , V = GU + gV 
or, in matrix notation, as 



Q 



U(x) 
V(x) 



Qc 



r Uo 

V 



N 



F(x) f(x) 
G(x) g(x) 



(15) 
(16) 

(17) 



Since D is constant, the determinant of the square matrix iV is a constant; its value, 
found by setting x = 0, is det iV = Fg — fG =1. It is usually more convenient to 
express Uo and Vo as a function of U(x) and V(x). Solving for Uq and Vo we obtain 
Q = MQ, where 

-/(*) 
F(x) 



M 



-G(x 

This matrix M is unimodular, |M 
Eqs. P) and (PUJ) as 

1 



1. Now we can find the characteristic matrix from 



M(x) 



cos 9 



cos (9 + Ex cos I 
—% sin (Ex cos 9 



—i sin (Ex cos 9) 
cos (9 — Ex cos 9) 



(19) 
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Figure 2. (a) Two opposite magnetic <5-function barriers, indicated by arrows (top 
figure) the vector potential showed by shaded (green) area, (b) The angles related to 
the propagation of an electron through this system (is shown in the bottom figure) . (c) 
Schematics of a periodic vector potential (shaded areas) and corresponding magnetic 
field indicated by the black arrows, (d) Arrangement of magnetic (5-functions that 
leads to a periodic area of vector potential with alternating sign. 
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Figure 3. The effective potential for k y < —P/2 and k y > P/2 for two different cases, 
(a) vector potential with non-zero average corresponding to Fig. Ufa), and (b) vector 
potential with zero average corresponding to Fig. [2tc). (c) Minimum (blue full curve) 
and maximum (red dashed curve) of the effective potential vs k y corresponding to the 
situation depicted in (a) for P = 0.1. 
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2.1. Bound states 

Regards to the average of vector potential we shall consider two different systems: one 
with zero average and the other with non zero average along the x-direction. First let 
us consider the magnetic field profile as shown in Fig. [2](a) for which the corresponding 
vector potential is 

A(x) =PB(x)Q(L-x), (20) 

where Q(x) = 0(x < 0), l(x > 0) is the theta function. This vector potential has a 
non- zero average, and the corresponding effective potential becomes (see Fig. [3(a)) as 

f k 2 y x<0 
V eff (k y ,x) ~ | (k y + P) 2 0<x<L . (21) 
{ k 2 x > L 

Here L is measured in the unit of magnetic length lg. There are two different cases 
which we have to consider. 

Case 1 for k y < —P/2: as shown in Fig. [3]^a) by the full red curve, we have a ID 
symmetric quantum well which, as is well-known, has at least one bound state (see also 
Fig. E(c)). For E 2 < k 2 the particle will be bound while for E 2 > k 2 we have scattered 
states, or equivalently the electron tunnels through the magnetic barriers. 

Case 2 for k y > —P/2: as shown in Fig. [3](a) by the dotted blue curve, the effective 
potential is like a barrier. We have a pure tunneling problem. With reference to Fig. [21 
Xx = and x<i = L, the solutions are as follows. For x < the wave function is 

il>(x) = Cer( 1 _ £ ), (22) 
where k = i?cosh^ and k y = i?sinh£, while for x > L it is 

il>(x) = De- M *( \ ). (23) 
In the middle region, < x < L, the wave function is given by 

^) = Fe*«( j, )+Qe-^(J_ w ). (24) 

With k! = [E 2 — (k y + P) 2 ] 1 / 2 = E cos 6. Matching the wave functions at x = and 
x = L leads to a system of four equations relating the coefficients C, D, F, and Q. 
Setting the determinant of these coefficients equal to zero, we obtain the transcendental 
equation, the solution of it gives the energy spectrum 

cos 9 cosh £ cos k'L + sin 9 sinh £ sin k'L = 0. (25) 

For the special value of k y = —P, sin 9 = 0, and we can rewrite Eq. ( 1251) as cos (EL) = 

or equally E n = + - J — . The resulting bound states, as a function of k y , are shown 

by the red full curves in Fig. H](a). The area of existence of bound states is delimited 
by the lines E = ~k y and E = —(k y + P). The number of bound states increases with 
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Figure 4. (a) Energy spectrum (red full curves on white background) and contour 
plot of the transmission (color background) through two magnetic, but with opposite 
direction <5- function barriers for L = 200 and P = 0.1. (b) the same as (a) but now 
for the configuration shown in Fig. [21(b) for a = b = c = 200 and P = 0.1. (c) 



Transmission vs energy through the magnetic <5-function barriers of in (a) k y 



-0.05. 



(d) same as (c) but now for the configuration shown in Fig. [2[b) for k y — 0. (e) 
Transmission vs k y through the system described in (a) for fixed E = 0.1. The red 
arrow lines indicate the position of the bound states, (f) same as (e) but now for the 
configuration shown in Fig. [D(b) for E = 0.15 



\k y \ which is also clear from the behavior of the min. and max. of the effective potential 
(see Fig. E^c)). No bound states are found for k y > —P/2 as is also apparent from Fig. 
[3](c) . For k y — > —P/2 the potential is shallow and only one bound state exists. The 
average velocity v n {k y ) along the y direction is given by 

/+oo 
dxj y (x), (26) 
-oo 

where j y = —i(U*V — V*U). From Fig. Hl^a) it is clear that these bound states move 
along the ^/-direction, i.e. along the magnetic barriers. Their velocity v y > —Vp is 
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negative for k y < —P but as the electron is approaching k y — > —P we have v y — > 0. 
For k y > —P the velocity v y > Vp is positive. This can be understood from the 
maximum and minimum of the effective potential which is shown in Fig. [3(c). The 
energy bound states can only exist between these two lines. Notice that the slope of 
min V e ff is negative for k y < —P while it turns positive for k y > —P which explains 
the k y dependence of the velocity. From Fig. H](a) it is clear there are two different 
classes of bound states. The bound state which follows very closely the k y = —P curve 
and extends to the region — P < k y < with energy close to zero has a wavefunction 
that is concentrated around the position of the two magnetic delta function and decays 
exponentially in the region < x < L. The wave function of the other bound states 
are concentrated in a region between the two magnetic delta-functions (i.e. like in a 
standing wave fashion) and decays exponentially outside this region. 

Next, we consider a structure with zero- average vector potential as shown in Fig. 
[2(b), with corresponding effective potential shown in Fig. [3(b). The effective potential 
for k y < —P/2 and k y > P/2, consist of a potential well and a potential barrier and 
therefore has at least one bound state. Thus we expect bound states for all k y with 
energy between E = -(+)(k y + P) and E = -(+)k y when k y < -P/2 (k y > P/2). The 
dispersion relation for those bound states are the solution of 

M 21 - Mi 2 - iM 22 e~t - iM n e f = 0, (27) 

where M is the transfer matrix for the unit shown in Fig. 12(d). These bound states are 
shown by the red full curves in Fig. E(b). Because of the spatial inversion symmetry 
of the vector potential the spectrum has the symmetry E(—k y ) = E(k y ). Notice that 
for — P < k y < P the lowest bound state has energy E m 0. For —P/2 < k y < P/2 we 
have two potential barriers and therefore no bound states. 



2.2. Reflection and transmission coefficients 

Consider a plane wave incident upon a system of two 5-function magnetic barriers, 
identical in height but opposite in direction, placed at x = and shown 
schematically in Fig. [2(a)- In this case the vector potential is constant for < x < L, 
zero outside this region, and homogeneous in the y direction. Below we derive 
expressions for the amplitudes and intensities of the reflected and transmitted waves. 

Let A, R, and T denote the amplitudes of the incident, reflected, and transmitted 
waves, respectively. Further, let #o be the angle of incidence and exit as shown in Fig. 
[2(b). The boundary conditions give 

U = A + R, U(L)=Te ikL , 
V = Ae ie ° - Re- i6 ° } V(L) = e i8 °e ikh T. 1 ' 



The four quantities Uo,Vq,U, and V given by Eqs. ( 1281) are connected by the basic 
relation Q = MQ; hence, with J = m' u + m' 12 e t9 ° and K = m' 21 + m 22 e t9 °, we have 

A + R=JT<**, 

Ae i9 ° - Re~ ido = KTe ikL , V ; 
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where m' • are the elements of the characteristic matrix of the medium, evaluated at 
x = L. From Eq. fl29|) we obtain the reflection and transmission amplitudes 

_ R _ Je ld ° - K _Z- 2e ~ ikhcos °o 

r ~A~ Je- ifl o + AT' * ~ A ~ Je-M* + K ' ^ ' 

In terms of r and t the reflectivity and transmissivity are 

ft=|rf, T=|t| 2 . (31) 

The characteristic matrix for a homogeneous vector potential is given by Eq. (fT9]) . 
Labeling with subscripts 1, 2, and 3 quantities which refer to the regions, respectively, 
I, II, and III of Fig. 12(a), and by L = x 2 — x\ distance between the magnetic ^-functions, 
we have {(3 = EL cos 9/) 

m' u = cos (9i + (5) / cos 9i, m' 22 = cos (9i — (3) / cos 9i, 

777,12 = — isin/3/ cos 9i, m 21 = — ism (3/ cos 9{. (32) 

The reflection and transmission amplitudes r and t are obtained by substituting these 
expressions in those for J and K that appear in Eq. (|30|) . The resulting formula can be 
expressed in terms of the amplitudes ri 2 ,ti2 and r 2 3,t 23 associated with the reflection 
at and transmission through the first and second "interface", respectively. We have 
e i6 ° - e^ 2cosfl 

ri2 e -i0 + e i0i ' tl2 e -i8 _|_ e ie s ' 

and similar expressions for r 2 3 and t 2 3. In terms of these expressions r and t become 
ri2 + r 23 e 2 ^ t 12 t 23 e^ 



1 + ri 2 r 23 e 2i/3 ' 1 + ri 2 r 23 e 2i/3 

The amplitude t of the transmission through the system is given by [H [9], [T2J [21] , 

2e~ ikL cos 9 cos 9i 

t = e-^cos^o + 9~) + 1] + e^[cos{9 - 9~) - 1] ' ( ' 3 '' j) 
where fc y = Esm.9$, and k y + P = Esin9i. This equation remains invariant under the 
changes E — ► — , ~~ ^ —9 ,9i — > — A contour plot of the transmission is shown 
in Fig. H](a) and slices for constant k y and E are shown respectively in Fig. S](c) and 
Fig. Hl(d). By imposing the condition that the wave number k x be real for incident and 
transmitted waves, we find that the angles 9q and 9i are related by 

$\n9 Q + P/E = sm9i. (36) 

Equation fl36l) expresses the angular confinement of the transmission elaborated in Refs. 
[H], [12], [2H], and [9]. Notice its formal similarity with Snell's law. Using Eq. ( 1361) we 
obtain the range of incidence angles for which transmission through the first magnetic 
barrier is possible 

-1 - P/E < sin# < 1 - P/E. (37) 

For the special value of the energy E = P/2 and #o i n the range —tt/2 < 9q < tt/2, we 
have 9i = 7r/2 while for E = —P/2 the result is 9{ = —n/2. Alternatively, we can put 
9i = ±7r/2 in Eq. ( |36l) and obtain, for P > 0, the result 

sin0Q = ±1 — P/E, (38) 
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where the +(— ) sign corresponds to E > {E < 0). A contour plot of the transmission 
as function of E and k y , obtained from Eq. fl35l) . is shown in Fig. H](a). In Fig. H^a) 
we distinguish three different regions. In the region between E = — (k y + P/2) and 
E = —k y , the wave vector of the incident wave is imaginary and they are evanescent 
waves. In this region k' is real and it is possible to find localized states. The k and k' 
for the second region between E = k y + P/2 and E = —(k y — P/2) are real and the 
electron can tunnel through the magnetic 5-barriers. In the blue shadow region between 
E = k y + P/2 and E = k y — P/2, k is real but k' is imaginary and solutions inside 
the barrier are evanescent and there is very little tunneling which becomes very quickly 
zero. The transmission probability \T\ = t ■ t* is equal to 1 for cos {2(3) = 1. In this case 
the energy becomes 



E„ = ± 



n 2 ix 2 /L 2 + (k y + P) 2 \ 1/2 , n = l,2,.... (39) 

The condition cos (2/5) = 1, or equivalently (3 = mr = Ehcos9 2 with n an integer, 
should be combined with that for the transmission to occur in the region delimited by 
the curves E = ±(k y + P) and E = ±k y . For example, in Fig. IU(a) for k y = —0.05 and 
< E < 0.2 we have 12 maxima. It is readily seen that with these parameters in Eq. 
(|39|) we find 12 different energies as shown in Fig. H{c). 

Fig. H](b) shows a contour plot of the transmission for the structure shown in Fig. EJ^d), 
which is symmetric around k y = 0. Notice that the number of resonances has increased 
substantially as compared to previous case which is due to the fact that we have twice 
as many magnetic barriers in our systems. 



3. A series of units with magnetic 5-function barriers 



3.1. N units 

We consider a system of N units, such as those shown in Fig. [21(a) and Fig. EJ^d) with 
periods L = a + b and L = a + b + c + d, respectively. The corresponding periodic 
vector potential is A(x) = A(x + nL) and the magnetic field B = B(x + nL), with 
n = 1, 2, N. The characteristic matrix for one period M(L) is denoted by 



M(L) 



m n m 12 

777.21 m 22 

On account of the periodicity we have 

M(NL) = M(L) ■ M(L)...M(L) 



(M(L)) 



JV 



(40) 



(41) 



N factors 

To evaluate the elements of M(iVX) we use a result from the theory of matrices, 
according to which the iVth power of a unimodular matrix M(L) is (wjv(x) = un) 



[M(L)\ 



N 



muUN-i - U N - 2 



m 12 u N -i 



(42) 
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with x = - Tr M and the Chebyshev polynomials of the second kind: 

^( X ) = sin[(iV + l)C]/sinC, (43) 

where 

C = cos" 1 x, (44) 

Here ( is the Block phase of the periodic system [33], which is related to the 
eigenfunctions of M. In the limit case of N — > oo, we have total reflection when ( 
is outside the range (—1, 1). 




Figure 5. Contour plot of the transmission (a) and Bloch phase (c) through N = 10 
magnetic ^-function barriers with a = 10, b = 10, and P = 1. (b) and (d) The same 
as in (a) and (c) for a = 5, b = 5, c = 5, d — 5, and P = 1, single unit. 




Figure 6. (a) and (b) Transmission vs energy through N — 1,5, 10 magnetic units 
of (5-function barriers shown on the left. The upper unit has a — 10, b — 10, P = 1, 
E = 1.5 and the bottom one a = b = c = d = 5, and P — 1 and E — 2.5. 
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3.2. Superlattice 

Here we consider a finite number N of lattice unit shown in Fig. [2(c). We set 

02 = Eb cos02, Pi = Ea cos 9\, P2 = 1/ cos #2, 



Pi = l/cos0i, h = a + b, X^ = 9 n ±f3 n . ^ 

The characteristic matrix M 2 (L) for one period is readily obtained, in terms of these 
quantities, as in Sec. II, and from that the characteristic matrix 'N\.2n{NL) of the 
multilayer system according to Eq. (|4ip . Its elements are 




it 

0.9-1 

01 




Figure 7. (a) The first two energy bands and contour plot of second (b) and first (c) 
band for a magnetic SL of 5-function barriers with a = 4, b = 4, and P = 1. 



Mn = s[cos A^ cos Xf — sin (3 2 sin Pi]un-i — un-2, 
M12 = — zsfcos sin/3i + sin/32 cos ]uat_i, 
M 2 i = — is [cos Xf sin/3 2 + sin/3i cosAj]itjv-i, 

M 2 2 = s[cos A^ cos Ajf — sin /3 2 sin /3i]iiiv-i — v>n-2, (46) 
where s = P2P1, ""at = un(x), and 

/ k 2 + k 2 + P 2 \ 

X = cos (fcia) cos (k 2 b) — — — sin (kia) sin (A; 2 &). (47) 

\ 2A; 2 fci J 

The reflection and transmission coefficients of the multi-unit system are immediately 
obtained by substituting these expressions into Eq. (13"U1) . The numerical results are 
shown in Figs. [51 [6] for finite superlattice with N = 10 units. Two different type of 
structures are considered as shown in the insets to Figs. El The transmission doesn't 
have k y — > —k y symmetry for the periodic system with magnetic delta up-down as is 
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M*L> k>L) k,(TiL) 

Figure 8. Dispersion relation (E vs k) for a standard electron in (a) and a Dirac 
electron in (b). The fixed values of k y are shown in the panels and L = 8 and P = 1 
(the energy for standard electron is measured in units of Tko c with uj c = y/ bBq / mc and 
all distances in Ib = \J ch/eBo). 




Figure 9. First energy band for (a) standard and (c) Dirac electron in SL of magnetic 
(5-function barriers with a = 100, b = 100, and P — 0.1. (b) and (d) corresponding 
contour plots of (a) and (c) (the energy for standard electron is measured in units of 
hu) c with u> c — \J eBo/mc and all distances in Ib — \J cfr/eBo). 

apparent from Fig. E^a). We contrast these results with the case in which we used 
an arrangement of magnetic cfe/ta-function as in previous structure plus another unit 
with opposite direction of magnetic delta function. As is clearly shown in Fig. [6(b), we 



Magnetic Kronig-Penney model for Dirac electrons in single-layer graphene 



15 




Figure 10. Dispersion relation for a (a) standard electron and a (b) Dirac electron. 
The period is L = a + b = 8 and the shaded (in green) regions are the lowest six allowed 
bands. The solid curves in both panels, the dash-dotted curves in (a) and the dashed 
ones in (b) show bound states for free electron (the energy for standard electron is 
measured in units of hu c with lo c — y/ cBq / mc and all distances in Is = y ch/eBo). 




12 14 5 6 7 8 12345678 
L L 



Figure 11. Energy vs period L — a + b for (a) standard electron and (b) Dirac electron 
with fixed k y — 4 (the energy for standard electron is measured in units of fko c with 
uj c = \J eBo/mc and all distances in Ib = \J 'ch/eBo). 



have k y — > —k y symmetry for the transmission probability through this structure. The 
transmission resonances are more pronounced, i.e., the dips become deeper, when the 
number of barriers increases for both types of units. But the gaps occur when the wave 
is mostly reflected. The position of these gaps, which are especially pronounced as N 
increases, can also be found from the structure of the Bloch phase £> as shown in Figs. 
|5|c) and (d). In Fig. [5] the bound states are shown by the blue solid curves that are 
situated in the area —k y —P/2 < E < —k y +P/2 in case (a) and in —k y —P/2 < E < —k y 
for case (b) plus an area located symmetric with respect to k y . Notice in Fig. [5](a) that 
several bound states merge into a resonant states at E = —k y + P/2. This is different 
from Fig. H^a) where each bound state becomes a resonant state at E = — k y . 
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Figure 12. (a) and (b): Lowest-energy band for a standard electrons and drift velocity 
vs k y for three different values of L. (c) and (d): First (black curves) and second (red 
curves) band for Dirac electrons and the drift velocity vs k y for three different values 
of L (the energy for standard electron is measured in units of hcu c with uj c = y eBo/mc 
and all distances in lg = \J ch/eB^). 



3. 3. Spectrum of a superlattice 

Lets take N — > oo. We can find the energy-momentum relation from the previous 
standard calculation [321 [33] by using 

cos{k x L) = ^TrM = X , (48) 

where M is the characteristic matrix of one period, which results into 

/fc 2 + k 2 + P 2 \ 

cos (kia) cos (k 2 b) — — — — !- ) sin (kia) sin (k 2 b) = cos k(b + a). (49) 

V 2 «2«i / 

With reference to the regions I and II shown in Fig. [2(a), we write k\ = [E 2 — k 2 ] 1 / 2 
and k 2 = [E 2 - (k y + P) 2 ] 1 / 2 an d s how the solution for E 2 > (k y + P) 2 in Fig. M 
Differences of Eq. ( l49l) from the corresponding result of Ref. [27], Eq. (JJj), for the 
case of Schrodinger electrons, is the term P 2 in the prefactor of the second term on the 
right-hand side and the linear E vs k spectrum instead of the quadratic one in Ref. [27J . 
If P is large the differences become more pronounced. Our numerical results for the 
energy spectrum are shown in Figs. [TJ [H [HJ CHI and [Til The results for standard and 
Dirac electrons show similarities but also important differences. The first band shows 
a qualitative difference near k y ~ 0, see Fig. [91 As Figs. [71 [9](a) and [9(b) show, the 
band behavior in the k = k x direction for fixed k y is constant and almost symmetric 
about k y = 0; the motion becomes nearly ID for relatively large k y . From the contour 
plots of Figs. Mjo) and (d), as well as from Fig. [7(c), we infer a collimation along the k y 
direction, i.e., v y oc dE/dk y w vf and v x ~ 0, which is similar to that found for a SL of 
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electric potential barriers [IB] for some specific values of the barrier heights. Also there 
are no gaps for k y « in Fig. [TOT a) but there are for the case of Dirac electrons as seen 
in Fig. [TOTb). This difference can be traced back to the presence of P 2 in the dispersion 
relation Eq. (H9l) when compare to the same equation for the standard electron. The 
even-number energy bands in Fig. fTTT b) are wider than those in Fig. [TIT a) and, as a 
function of the period, the energy decreases faster for Dirac electrons. This behavior of 
the bands for Dirac electrons is very similar to that for the frequency u vs k y or L in 
media with a periodically varying refractive index [29] . This is clearly a consequence 
of the linear E — k relation. Notice the differences between the lowest bands shown in 
panels (a) and (b) in Fig. [12] and in particular the difference between the corresponding 
drift velocities as functions of k v . 



4. A series of 5-function vector potentials 

In the limit that the distance between the opposite directed magnetic barriers decrease to 
zero the vector potential approaches a 5-function [24J . We consider a series of magnetic 
5-function vector potentials A(x) = Y^L-oa Aq6(x — nL) as shown in Fig. 11(a). First 
we consider a single such potential, that is zero everywhere except at x = 0. We start 
with Eqs. ([6]) and ([7]) which becomes now 

- i [d/dx + k y + £ B A S(x)] (fi b = E(p a , (50) 
— i [d/dx — k y — £ B A 5(x)} (fi a = E(p b . (51) 

The solutions are readily obtained in the form 

| A cos (ex cos 9) + B sin (ex cos 9) , x < 0, 

\ C cos (ex cos 9) + D sin (ex cos 9) , x > 0, 

| — i {B cos (9 + ex cos 9) — A sin (9 + ex cos 9)} , x < 0, . 

\ — z {.Dcos (9 + ex cos 9) — C sin (9 + ex cos 9)} , x > 0. 

Integrating Eqs. (1501 . (!5Tj) around gives 

/•o+ <-o+ 

— i / [d/dx + (A™, + £bAo5(x))} fadx = E / 4> a dx, (54) 

jo- jo- 

—i / [d/dx — (k y + £bAq5(x))] 4> a dx =E <pbdx (55) 
Jo- Jo~ 

and 

<M0") = # fe (0 + ), a (O + ) = #a(0"). (56) 

We now consider the entire series of 5-function vector potentials shown in Fig. fTST a) and 
use Eqs. ( )54l) and the periodic boundary condition ^/(O) = e lkxL ^ n(L). The resulting 
dispersion relation for the superlattice is 

cos (kL) = 277(1 + f? 2 )" 1 cos(k x L), (57) 
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Figure 13. (a) A series of ^-function vector potentials, (b) Dispersion relation for the 
system shown in (a). 



where n — 1 + (-bAq. From Eq. ( I5"7j) we can find the energy spectrum as 

E n (k x , ky) = it 



2nn + k?.+ \— cos" 



2?/ 



1+77' 



cos (k x L) 



We can define 



(1/L) cos^ 1 f2?7(l + r/ 2 )" 1 cos{k x L) 



and obtain 



£ = ±[2rar + k 2 + s 



211/2 



(58) 



(59) 



(60) 



The energy bands around the Dirac point are plotted in Fig. [T2Tb). Notice that: 1) 
there is an opening of a gap at the Dirac point, 2) the motion is strongly ID, i.e. along 
the ^-direction, and 3) higher subbands have a smaller dispersion. 



5. Concluding remarks 

We developed a magnetic Kronig-Penney model for Dirac electrons in graphene. The 
model is essentially a series of very high and very narrow magnetic ^-function barriers 
alternating in signs. The treatment of the transmission through such a series of barriers 
followed closely the one developed in optics for media in which the refractive index varies 
in space [2S1 ES]- We contrasted a few of the results with those for standard electrons 
described by the Schrodinger equation [27] . 

In several cases the energy spectrum or the dispersion relation were obtained 
analytically, cf. Eqs. (25), (39), (47), (48), and (57), largely due to the simplicity 
of the model and the adapted method from optics. For only two magnetic 5-function 
barriers, opposite in sign, we saw several bound states, whose number increases with 
\k y \, and a reduction of the wavevector range for which tunneling is possible, cf. Fig. 
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4(a). This is in line with that reported earlier for single [7j and multiple [H] barriers. 
The reduction becomes stronger as we increase the number of barriers, cf. Fig. 4(b). 
We also made contact with Snell's law in optics, cf. Eq. (36): the term P/E represents 
the deviation from this law. 

An important feature of the superlattice results is a collimation of an incident 
electron beam normal to the superlattice direction at least for large wave vectors. As 
easily seen from Figs. 7 and 8, for \k y \ > 2 we have v x oc dE/dk x « for the first 
three minibands in the middle panels and nearly five minibands in the right panels. 
This occurs for both standard electrons and Dirac electrons but notice an important 
difference for \k y \ « shown clearly in Fig. 9. This collimation is similar to that 
reported in Ref. [16] for superlattices involving only electric barriers but with somewhat 
unrealistic large barrier heights. 

It is also worth emphasizing the differences and similarities in the first two 
minibands and the corresponding drift velocities as functions of k y for different periods 
L and constant k x as shown in Fig. 12. Notice in particular the resemblance between the 
drift velocities in the lowest miniband for standard electrons and the second miniband 
for Dirac electrons. 

Given that ferromagnetic strips were successfully deposited on top of a 2DEG in 
a semiconductor heterostructure [23], we hope they will be deposited on graphene too 
and that the results of this paper will be tested in a near future. 
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